Vol.11, No.1, Spring 2022, pp.163-183 ISSN (Print):2322-1682 ISSN (Online):2383-3076 


#یرانی ژد 
Se,‏ یرانی oe‏ 


Journal of Geography 


and Environmental Hazards Girth > k 
دانشگاه فردوسی مشهد‎ Serer een 


Ferdowsi University of Mashhad 


RESEARCH ARTICLE DOI: 10.22067/geoeh.2021.71958.1096 Open access 


Comparative Study of the Best Methods for Estimating Maximum 
Annual Precipitation in Selected Stations of Haraz River 


© 
Abbas Gholami“ 


Assistant Professor, Faculty of Engineering ,Shomal University,Amol, Mazandaran, Iran 


Received: 31 August 2021 Revised: 8 October 2021 Accepted: 21 October 2021 


Abstract 


The purpose of this study was to compare the statistical parameters of maximum annual 
precipitation in the watershed of Haraz River located in the central part of the Alborz 
Mountains. After reviewing the data of precipitation in Haraz watershed and considering lack of 
long-term statistics in the study area, Panjab-Haraz, Karsang-Haraz and Mahmoudabad-Haraz 
stations in different parts of the river from mountain to plain in 44-year statistical period from 
the water year of 1968-69 to 2011-2012 were selected. The annual precipitation frequency 
analysis of the selected stations was performed and the results of the Residual sum of squares 
(RSS) test were extracted by the method of ordinary moments and the maximum likelihood at 
this stage to be compared with the linear moments method. Residual sum of the squares in the 
linear moment method were also calculated by a computer program for the studied distributions. 
The results showed that in general, normal moment and maximum likelihood methods in 
selected stations have better fit and less RSS. The results of scoring and ranking the best 
selected method and distribution also showed that in all three selected stations, the maximum 
likelihood method with the minimum RSS was the best fit with the annual maximum 
precipitation data. 
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جغرافیا و مخاطرات محیطی. سال یازدهم. sled‏ چهل و یکم بهار ۶۰۱ صص ۱۳۳-۱۸۳ 


مقاله پژوهشی 
بررسی مقایسه‌ای بهترین روش برآورد بارش حداکثر سالانه در ایستگاه‌های منتخب رودخانة هراز 
® عباس غلامی "- استادیار دانشکده فنی و مهندسی,دانشگاه شمال.آملمازندرانایران 
تاریخ دریافت: ۱۰۰/۷۹ تاریخ بازنگری: ۱۶۱۰/۷/۱۱ تاریخ تصویب: ۱۶۱۰/۸/۲۹ 


حکیده 


با توجه به اهداف پژوهش حاضر پس از بررسی آمار و اطلاعات بارش‌های حوزه آبخیز هراز و با توجه به 
کمبود آمار طولانی‌مدت در منطقه موردمطالعه» درنهایت ۳ ایستگاه پنجاب- هرا کرسسنگ- هراز و 
محمودآباد- هراز در بخش‌های مختلف این رودخانه از کوهستانی تا جلگه‌ای در دوره مشترک آماری ٤٤‏ 
تساه اسان ای ۱۳۱۳۵۸ ۱۳۹۱-۵۱ تشاب ts‏ سس از Litas‏ کیرات gla fol‏ سالا 
ایستگاه‌های انتخاب شده نتایج آزمون حداقل مربعات به روش گشتاورهای معمولی و بیشینۀ درستنمائی 
برای مقایسه با روش گشتاورهای خطی استخراج شدند. مجموع مربعات باقیمانده در روش گشتاورهای 
خطی نیز با برنامه کامپیوتری نوشته شده برای توزیع‌های موردپژوهش, محاسبه گردیدند. با توجه به نتایج 
این پژوهش و با استفاده از روش نمره دهی و رتبه‌بندی مشخص می‌شود در ایستگاه پنجاب- هراز روش 
بیشینۀ درستنمائی دارای کمترین مجموع مربعات خطاء یعنی مجموع نمره 1 دارای بهترین برازش با 
داده‌های بارش سالیانه بوده است و روش گشتاورهای معمولی و خطی با دارا بودن iya‏ کل ۱۵ به‌صورت 
مساوی در رتبه‌های بعدی قرار دارند. در ایستگاه کرسنگ- هراز روش dade‏ درستنمائی دارای مجموع 
نمر ۱۰ و بهترین برازش با داده‌های بارش سالانه بوده است و پس‌ازآن به ترتیب روش گشتاور معمولی با 
مجموع VE Spa)‏ و روش گشتاور خطی با مجموع نمرةٌ ۱۵ قرار دارند. در ایستگاه محمودآباد هم روش 
بیشینۀ درستنمائی با به دست آوردن مجموع نمرة ۳ بهترین روش منتخب در این ایستگاه بوده است و 


پسازآن روش گشتاور معمولی و روش گشتاور خطی به ترتیب با مجموع نمرات ٩‏ و ۱۵ قرار دارند. 
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Lal ae ۳4‏ و مخاطرات محبطی شمارة اول 


۱- مقدمه 

حوادث بارندگی و در نتیجه آن سیل‌ها می‌توانند سبب خسارت‌های فراوانی در نقاط مختلف دنیا شوند. مدیریت 
دشت‌های سیلابی و طراحی برای فعالیت‌های JES‏ سیلاب. سدهاء پل‌ها و دیگر پژوهش‌های موردنیاز در امکان 
وقوع يا احتمال هر یک از حوادث منعکس می‌شود. (استدینگر" و همکاران» ۱۹۹۳). مشکل عمده تجزیه‌وتحلیل 
نوسانات بارندگی کمبود آمار و داده‌های هیدرولوژی در bE‏ مختلف است. در عمل توزیع احتمال واقعی داده‌ها در 
یک مکان. مجهول است و با روش پارامتری ساختار تابع چگالی احتمال معلوم می‌شود و فقط باید پارامترها را با 
استفاده از داده‌ها برآورد نمود. برای تخمین پارامترها از روش‌های مختلفی استفاده می‌شود. (فیضی و اسلامیان" ۱۳۸۶). 
تجربیات طولانی نشان داده است و بر اهل‌فن پوشیده نیست که بدون استفاده از قوانین آمار و احتمالات و قبل از 
تجزیه‌وتحلیل داده‌های هیدرولوژی و هواشناسی استفاده از داده‌های ple‏ در طراحی ساختمان‌های آبی خالی از 
اشکال نیست و ممکن است خطرناک ب‌اشد. (سیمافر " ۱۳۷۰). یکی از روش‌های آماری استفاده از انواع توزیع‌های 
فراوانی یا آماری جهت مشسخص کردن فراوانی وقوع یا تعداد دفعاتی است که به‌طور متوسط یک متغیر در مدت‌زمان 
مشخص اتفاق می‌افتد(علیزاده » ۱۳۷۲). با در نظر گرفتن اینکه دوره‌های آمارمرداری کوتاه می‌باشند سعی می‌گردد 
تا منحنی به دست آمده را با یکی از توزیع‌های آماری تطابق داده و بهترین منحنی را انتخاب کنند. سپس با مشخص 
بودن منحنی تطابق یافته» دبی سیلاب‌ها یا کم‌آبی‌ها با دوره‌های برگشت متفاوت تعبین می‌گردند. تطابق توزیع 
فراوانی داده‌های مسوجود با یکی از توزیع‌های آماری این حسن را دارد که آمار محدود موجود را می‌توان توسعه داد 
و احتمالات بسیار کم را نیز از روی امتداد منحنی‌های توزیع تلوری به دست آورد و ضمناً می‌توان توزیع حاصله را با 
دو و یا سه پارامتر بیان نمود (مهدوی؟ ۱۳۷۱). در این پژوهش برای برآورد پارامترهای آماری از سه روش گشتاور 
معمولی بیشینه درستنمائی و روش جدیدتر روش گشتاور[ (گشتاورهای خطی) استفاده شده است. لازم به توضیح 
است که پارامترهای آماری محاسبه شده از روش گشتاور 1 نسز می‌توانند مانند روش گشتاور خطی و بیشینة 
درستنمائی برای محاسبة بارش‌های حداکثر سالانه مورداستفاده قرار گیرند و درنهایت با استفاده از آزمون‌های برازش 
نکوئی» بهترین توزیع برای ایستگاه‌های مختلف به دست می‌آید. در این پژوهش نیز کارایی این سه روش برای 
محاسبه پارامترهای آماری بارش حداکثر ساانه مقایسه شده است که این موارد ازجمله اهداف اصلی تحقیق 

می‌باشند. سؤال اصلی این پژوهش به‌صورت زیر است: 
Stedinger‏ 1 
Feizei and Islamyan‏ 2 
Simafar‏ 3 


4 Alizadeh 
5 Mahdavi 


سال یازدهم بررسی مقایسه‌ای بهترین روش برآورد بارش حداکثر .... me‏ 


استفاده از کدام روش محاسبه پارامترهای آماری و کدام توزیع در ایستگاه‌های انتخاب شده دارای حطای کمتر و 
دقت بالاتری در محاسبه بارش سالیانه این منطقه است؟ 

پاسخ صحیح به این سژال ازجمله Glial‏ اصلی این تحقیق است. 

فیضی و اسلامیان (۱۳۸۶) تحلیل فراوانی منطقه‌ای حداکثر بارش ماهانه را با. استفاده از روش گشتاور خطی در 
۰ ایستگاه حوضه زاینده‌رود موردبررسی قرار دادند و مشاهده کردند که aly‏ توزیع مقادیر حدی تعمیم یافته و 
پیرسون نوع سه بهترین برازش را به داده‌های حداکثر بارش ماهانه دارند.حقیقت جو" (۱۳۸۱) روش‌های غیرپارامتری 
را برای تحلیل فراوانی پنج ایستگاه قدیمی اصفهان بوشهر جاسک» تهران و مشهد انجام داده و نتایج را با روش 
گشتاور معمولی مقایسه کرد. 

در ابتدا وقتی گشتاورهای L‏ (هاسکینگ 194( به‌وسیلهٌ محققان مختلف توصیه شد آن‌ها از آن برای برازش 
نکوئی توزیع‌های احتمالاتی مختلف در نمونه‌های منطقه‌ای جریان‌ها و بارش‌ها اسستفاده کردند. پیرسون * (1۹۹۵) 
نمودارهای گشستاور L‏ را برای نزدیک به ۵۰۰ سری جریان حداقل سالانه در نیوزیلند بکار برد. او پی برد که تنها 
مدل‌های ۲ پارامتری برای مشاهده روابط گشتاور 1 در نیوزیلند کافی نیست. با استفاده از اطلاعات جریانات متوسط 
سالانه در ۹۶۷ منطقه در اطراف جهان» مک ماهون* وهمکاران (۱۹۹۲) و فینلیسون و مک ماهون' (۱۹۹۲) ثابت 
کردند که تقریباً ۰ رودخانه‌های جهان توزیع‌هایی دارند که احتلاف معنی‌داری از توزیع نرمال ندارند (براس‌اس 
آزمون چولگی با حدود اعتماد ۵./). یک آزمایش محدودتر ding‏ مک ماهون" و همکاران (۱۹۹۲) نشان داد که تنها 
بالغ بر ۸۳۱ مناطق در استرالیا و آفریقای جنوبی تقریباً نرمال هسستند که با ۷۱/ مناطق در رابطه با داده‌های جهانی 
قابل‌مقایسه است. مشابه آن در مطالعات پیش‌تر و محدودتر با استفاده از ۱۳۷ و ۱۶۰ حوزه آبخیز در اطراف جهان 
کالینین ‏ (۱۹۷۱) و یوجویچ* OAW)‏ به ترتیب پی بردند که جریانات سالانه در بالغ بر ۰ مناطق تقریباً دارای توزیع 
نرمال هستند. مارکوویچ " OAIT)‏ از روش آماری برازش نکوئی حداقل مربعات برای برازش توزیع‌های LN 2 ٩‏ 
LN 3‏ ]0۸1۷ و3 ۴ در توالی جریانات متوسط ساانه در LEN‏ منطقه در غرب ایالات‌متحده استفاده کرد. مارکوویچ 
استفاده از توزیع گام‌ها را توصیه کرد اگرچه ننایج نشان داد که توزیع‌های 1N 3 5LN2 P3‏ ب‌صورت متغیرهای 


1 Feizei and Islamyan 

2 Haghighatjoo 

3 Hosking 

4 Pearson 

5 McMahon 

6 Finlaysonand Mc Mahoon 
7 McMahon 

8 Kalinin 

9 Yevjevich 

10 Markovich 


Lal ae wv‏ و مخاطرات محبطی شمارة اول 


گرفته شده به‌وسیله مارکوویچج است. لف و هاردیسون" (VAT)‏ از پلات‌های احتمالی برای حدس O25‏ برازش نکوئی 
توزیع‌های نرمال. لوگ les‏ ویبول برای جریانات متوسط سالانه در YY‏ منطقه منابع آب در ایالات‌متحده استفاده 
کردند آن‌ها توزیع نرمال را برای حوزه‌های آبخیز نیوانگند. اوهایو» تنسی. می‌سی‌سی‌پی. میسوری. کلرادو و شمال 
غرب اقبانوس آرام انتخاب کردند توزیع لوگ نرمال را برای حوزه‌های آبخیز دلاور» هادسون» جنوب شرق» ریوگرانده 
بالایی» پسوس و اقیانوس آرام مرکزی و توزیع ویبول را برای باقیمانده حوزه‌های آبخیز انتخاب کردند. ووگل و 
همکاران (1۹۹۵) با استفاده از نمودارهای گشتاور L‏ و هاسکینگ" (1۹۹۰) با آزمایش نرمال‌بودن گشتاور gL‏ 
as task ely, alas‏ قاب OLAS Soe ON‏ که Ul go‏ بر سظ ails‏ در ۱4۳ مط در شمان 
شرق ایالات‌متحده برازش نسبتاً عوبی با یک توزیع نرمال دارند. ایگلسون* (1۹۷۸) توزیع احتمالی جریان سانانه را با 
فرض یک تابع ساده بارش- بازده به‌وسیله رگبارهای رسیده بر طبق روش پوآسن و عمق رگبار مطابق با یک توزیع 
گاما نتیجه گرفت. (VW) Ke pb‏ یک مدل تعادلی آبی ساده بارندگی هرز آب را بکار برد که شامل ۲ معادله است: 
یکی برای جریان رودخانه‌ای و یکی برای ذخیره آب زیرزمینی. سالاس و اسمیت" (VAM)‏ خاطر نشان کردند که این 
مجموعه مدل بارندگی- رواناب شاید به‌مانند یک میانگین متحرک رگرسیونی تفسیر شود (مدل تجربی جریان 
رودخانه‌ای به همراه روش‌های معمول براساس سری‌های بارش). در یک ارزیابی اخیر ازسری‌های زمانی بارش 
متوسط سالات در کل ایالاتمتحده با استفاده از نمودارهای گشتاور D‏ گوتمن " و همکاران (1۹۹۳) مدل‌بندی pt‏ را به 
همراه یک توزیع PY‏ توصیه می‌کنند. مارکوویچ" )1900( توزیع گاما را برای جریان‌های سالانه توصیف می‌کند و 
همین‌طور cle‏ تعجب نیست که بیشتر سری‌های جریان سالانه که به‌وسیله کالینین "و همکاران (۱۹۷۱) یوجویج V‏ 
OAT)‏ و ووگل" و همکاران )480( ارزیابی شد برازش Gs‏ خوبی با یک توزیع نرمال دارند. بر Gb‏ 
پژوهش‌های لل " و همکاران ( 1۹7۳ با افزایش اندازه نمونه» انتخاب روش پارامتری مناسسب» آسان‌تر اسست و 
هم‌چنین» در این گونه نمونه‌ها پارامترهای برآورد شده توأم با خطای کم‌تری هسستند. در تحقیقی که بوستالا" و 


همکاران (۲۰۱۱) روی چهار زیر حوزه در ایرلند انجام دادند به بررسی SE‏ تغییر اقلیم روی دبی سیلاب پرداخته و 


1 Lof and Hardison 
2 Vogel 

3 Hosking 

4 Eagleson 

5 Fiering 

6 Salas and Smith 
7 Guttman 

8 Markovich 

9 Kalinin 

10 Yevjevich 

11 Vogel 

12 Lall 

13 Bastola 


سال یازدهم بررسی مقایسه‌ای بهترین روش برآورد بارش حداکثر .... YW‏ 


ابتدا با آنالیز حساسیت دبی‌های سیلابی نسبت به پارامترهای دماء بارش و تبخیر و تعرق دریافتند که دبی پیک بیشترین 
حساسیت را نسبت به بارش دارد. سپس با استفاده از گشتاورهای خطی به جداسازی مناطق همگن پرداخته و به تعیین 
دبی پیک در آینده با استفاده از پارامتر بارش پرداختند بتول ! (۲۰۱۰) در تحقیقی که در یکی از حوزه‌های آبخیز ترکیه 
انجام داد به بررسی اثر ایستگاه‌های غیرهمگن(ناجور) در یک منطقه و تأثیری که روی پیش‌بینی دبی پیک با دوره 
بازگشت‌های مختلف می گذارد پرداخت به این منظور ابتدا با استفاده از گشتاورهای خطی به جداسازی مناطق همگن 
پرداخته و دو شاخحص ناجوری را با استفاده از توزیع ٤‏ پارامتره US‏ به دست آوردند نتایج حطای مربعات استاندارد 


فراوانی سیل نشان داد که مناطق با شاخص همکنی قوی بسیار دقیق‌تر از مناطقی با شاخص همکنی ضعیف بودند. 


۲- مواد و روش‌ها 

۱1-۲- موقعیت جغرافیایی منطقة مورد مطالعه 

رودخانه هرازه رودی است در استان مازندران در شمال lel‏ که از دره لار در جنوب oS‏ دماوند سرچشمه گرفته 
و به دریای خزر می‌ریزد. پس از سرازیر شدن رود هراز از دره oU‏ رواناب‌های دره قاضی US‏ و جنگل‌دره از کوه‌های 
شمیم‌کوه» آب زیارو و لاسم از کوه‌های غربی قزقان‌چای فیروزکوه و آب رودهای جاری از دهکده‌های دلارستاق, نیز 
به این رود وارد می‌شسوند. آب زراعی کشاورزان آمل» فریدونکنار بخشی از بابل و نور نیز از این رودخانه تأمین 
می‌گردد. مهم‌ترین شسهرهایی که بر روی مخروط افکنه رودخانه هراز توسعه يافته‌اند عبارتند از آمل و فریدونکنار. 
شیب رودخانه هراز در محدوده کوهستانی بسیار متغیر است. شیب رودخانه هراز از مرز کوهستان تا شمال شهر آمل 
۳ در هزار و در محدوده شهر آمل ۷ در هزار است. رودخانه هراز در محدوده کوهستانی» در دره‌ای نسبتاً باریک 
جریان دارد. با نزدیک شدن به پهنه ساحلی» پهنای آن افزايش يافته و رودخانه به دلیل افت شیب بستر به چندین 
شاخه تفکیک می‌گردد مهم‌ترین سرچش مه رود She‏ از قله 1۳۷۵ متری Wy‏ گردن در VW‏ کیلومتری جنوب 
شرقی چالوس با نام لار جاری می‌شسود و پس از دریافت آب‌های مهم دیگری (مثل: آب سفید. الرم. آب چهل بره 
آب سیاه پلاس» آب امام پهنک» سه سنگ. دیو آسیاب ورارود دلی چای و چشمه ملک) و عبور از ده پلور با نام 
هراز به سمت دریای مازندران جاری می‌گردد. البته در راه باز رودها و چشمه‌هایی به آن می‌پیوندند که مهم‌ترین 
lol‏ لاسم. تلخ رود آب رزان» آب مشک انبار» پنجاب (نمارساق) هراز» آب پردمه شیر کلارود. چلورود و (Skee he‏ 
فسا طرل رود هراز teal els Wo‏ و از کار GS wiles ll alt gut‏ وباک از این شیر 
لاریجان نام دار که نام خود را از روی سرچشمه اصلی هرازه لار گرفته است. در مسیر این رودخانه دو سد زیبای 
لار و هراز قرار دارد. شسکل(۱) موقعیت ایستگاه‌های هیدرومتری موردبررسی در بخش مرکزی استان مازندران را 
نشان می‌دهد. 


1 Betül 


Lal ae ۱۷۹‏ و مخاطرات محبطی شمارة اول 


شکل ۱- موقعیت ایستگاه‌های هیدرومتری موردبررسی در بخش مرکزی استان مازندران 


۲-۲ روش انجام پژوهش 
۱-۲-۲- ایستگاه‌های هیدرومتری منطقه 


پس از بررسی آمار بارندگی کلیه ایستگاه‌های تحت پوشش حوزه آبخیز هراز از نظر مدت آماری: یکنواختی آمار 
سال‌های مختلف و عدم وجود سال‌های بدون آمار در منطقه, ۳ ایستگاه با مشخصات ذیل برای بررسی انتخاب 


شدند. جدول (۱) جزئیات ایستگاه‌های باران‌سنجی انتخاب شده در این پژوهش را نشان می‌دهد. 


حدول ۱- مشخصات ایستگاه‌های باران‌سنحی انتخاب شده در منطقه موردمطالعه 


رودخانه مشخصات جغرافیایی 
ایستگاه طول عرض ارتفاع 
پنجاب- هراز ay. ۳۹۲-۱۵۰ o-٤‏ 
کره سنگ-هراز ۵۲-۲۲-۰۵ ۳-۵ Yo‏ 
محمودآباد-هراز ۵۲-۱۵-۳۱ ۳-۷-۷ -Yo‏ 


با توجه به مطالعات انجام شده بر روی آمار ایستگاه‌های منتخب و با توجه به استقرار این ایستگاه‌ها در مکان‌های 
مختلف حوزه و با توجه به اینکه آزمون‌های مورداستفاده در این تحقیق نیاز به وجود یک سری آماری بارشی مداوم 
دارد» بنابراین سال آماری 44- ۱۳۶۸ تا سال آماری ۱۳۹۰-۹۱ به مدت 4۶ سال به‌عنوان یک دوره آماری مشسترک 


مختلف و تأثیر این موضوع بر نوسانات بارش ناشی از مناطق دارای برف باران پا برف و باران بود. با توجه به اینکه 


سال یازدهم بررسی مقایسه‌ای بهترین روش برآورد بارش حداکثر .... ۳1۹ 


توزیع بارندگی ماهانه فصلی و یا سالیانه در طول سال می‌تواند نقشی تعیین کننده در محاسبات هیدرولوژیکی و دبی 
فصلی و سالانه در جهان منجر به ظهور رژیم‌های ویژه بارشی شده است که می‌توان به رژیم‌های استوایی. گرمسیری» 
مدیترانه‌ای اشاره کرد. 

۲-۲-۲- روش‌های مختلف برآورد پارامترهای توابع توزیع آماری 

۱-۲-۲-۲- روش بيشينة درستنمای ) (Maximum Likelihood Method‏ 
بر اساس بیشینه درست- نمایی برای رخدادهای مشاهداتی به دست می‌آید پیفانزل (1۹۹۶). اصل iuiu‏ درستنمایی 
ماکزیمم می‌کند که برای یک توزیع یا تابع چگالی احتمال (....,8 pla,‏ که در آن » و 8 پارامترهای توزیع بوده و 
باید برآورد شوند احتمال به دست آوردن مقداری از × مانند ,× متناسب با (.... , 8 + phia‏ و احتمال توم 
L‏ برای به دست آوردن نمونه‌ای از 1 مقدار 260 ,..... ,262 ,261 که متناسب است با حاصل‌ضرب زیر: 

L=[[ p(x, a, Boan) (\) 
i=] 

ماکزیمم شود. این ماکزیمم از مشتقات نسبی مآنسبت به هر یک از پارامترها که مساوی صفر قرار گرفته باشند به 
دست می‌آیند. اغلب برای Sale‏ از Lab‏ بجای L‏ در محاسبات استفاده می‌شود. 
۲-۲-۲-۲- روش کمترین مربعات ) (M ni mumsquare net hod‏ 

روش برآورد کمترین مربعات عبارت از برازش یک تابع نظری با یک توزیع تجربی اسست. پس‌ازاین برازش 
مجموع مربعات تمام انحراف‌های تقاط مشاهده شده از تابع برازش را مینیمم می‌کنند. پس‌ازاین برازش تابعی 
به‌صورت: 


P= تاه‎ eae) (Y) 


S=) (yi =F Gis رو‎ ,..)) G) یا‎ 
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را Lb‏ مینیمم سازيم. که در آن و yi‏ مختصات نقاط مشاهده شده 0 و 8 پارامترها و N‏ حجم نمونه است. برای 


(a)‏ 0= ۱9۸ ده 


i=l 


از این مشتقات نسبی به تعداد پارامترهایی که باید برآورد شوند معادله به دست می‌آید. 


9 
(Graphical Method ) روش گرافیکی‎ -۳-۲-۲-۲ 


پارامتره روی منحنی m‏ نقطه اختیار می‌کنيم تا m‏ معادله برای حل کردن به دست آید. این فرایند را می‌توان با امستفاده 
از pl il‏ کاغذ گراف و استفاده از تغییر مختصات تا حدی که برازش یک ط راست امکان‌پذیر شود ساده‌تر کرد. 


(Ordinary Moments Method ( روش گشتاورهای معمولی‎ -۶-۲-۲-۲ 


قدیمی‌ترین و گسترده‌ترین روش شناخته شده برای برازش فرکانس توزیع‌های مشاهداتی. روش گشتاورهاست 
(ووگل و همکاران ۱۹۹۳). در روش گشتاورها Vales‏ کلی برای محاسبة گشتاور alt AS yo‏ حول مبداً یک توزیع» 
P (x)‏ به کار می‌رود: 


0 


4. = fx -p(2).dx ۹2 


00 


یا dole‏ متناظر برای گشتاورهای مرکزی توزیع مورداستفاده قرار می‌گیرد: 
u. = 7 y .p(x).dx (A)‏ 


که در آن ME‏ گشتاور مرتبة اول حول مبداً است. روش گشتاورهاء ilih‏ بین گشتاورها و پارامترهای توزیع را 


مشخص می‌سازد. عموماً روش گشتاورها بوسیلة ایجاد Wales‏ گشتاورهای تثوریکی قابل کاربرد است: 


E[x]‏ بر 
=Var|x|= El(x = uy]‏ ?5 
Alle]‏ 
?8 
x-£le-)‏ 
8l‏ 


سال یازدهم بررسی مقایسه‌ای بهترین روش برآورد بارش حداکثر .... wi‏ 


alsa (Vs) 8 =| 2S Gin) 0۱ 
2 2 
K= 1 1 n ۱ 4 1 1 11 ۱ 3 
7 Xhi m) (VY) G=— -X (xi-m) (۳) 


که !۸ , KV‏ به ترتیب میانگین» واریانس» چولگی و کشیدگی تعریف می‌شوند a > , G , 8” ms‏ 
ترتیب میانگین نمونه» واریانس» چولگی و کشیدگی تعریف می‌شود. این موارد جزء اصولی‌ترین مشسخصات آماری 
می‌باشند معمولاً تمایل مرکزی توسسط میانگین بیان می‌گردد و برای بیان پراکندگی و درجة تقارن از واریانس ‏ و 
ضریب چولگی " و ضریب کشیدگی " استفاده می‌شود. در غالب کارهای هیدرولوژیکی امکان دستیابی به مشخصات 
آماری کل جامعه وجود ندارد. به همین دلیل در بررسی این‌گونه جوامع» مشخصات آماری نمونه‌هایی از جامعه 
حساب شده و براساس Ol‏ مشخصات آماری کل جامعه تخمین زده می‌شود. هرچه بر تعداد اعضای جامعۀ نمونه 
افزوده شود مشخصات آماری حساب شده به مقدار واقعی OT‏ در کل جامعه نزدیک‌تر خواهد شد. تمایل مرکزی 
داده‌ها توسط سه کمیت آماری بیان می‌شود. متوسط حسابی داده‌ها یا میانگین قابل‌اعتمادترین شاخص pled‏ مرکزی 
آن‌هاست. این شاخص از daly‏ (۱۰) به دست می‌آید که در آن n‏ تعداد نمونه‌ها و ۳ میانگین داده‌هاست. هر گاه 1 به 
سمت بی‌نهایت میل LS‏ مقدار × نیز سمت میانگین جامعث ل میل خواهد کرد میزان پراکندگی داده‌ها توس ط 
واربانس آن‌ها مشخص می‌گردد این معیار معرف میانگین مربع انحرافات داده‌ها از میانگین اسست و از (VV) ib‏ 
تعیین می‌شود. جذر واریانس, انحراف معیار* نامیده می‌شود و در مسائل آماری مورد baled‏ فراوان قرار می‌گیرد. هر 
چه انحراف معیار داده‌ها کمتر باشد پراکندگی آن‌ها در اطراف میانگین کمتر خواهد بود. درحالی که ميانگین و 
واریانس. گشتاورهای اول و دوم داده‌ها را تشکیل می‌دهند» برای بررسی تقارن آن‌ها از گشتاور سوم استفاده می‌شود 
و ضریبی به نام ضریب چولگی تعریف می‌شود که آن‌هم از رابطۀ (۱۲) به دست می‌آید. اگر یک توزیع نسبت به یک 
توزیع متقارن انحراف داشته باشد آن توزیع را توزیع چوله می‌نامیم. گشتاور چهارم حول میانگین کشیدگی است که 
بلندی یا کوتاهی آن نسبت به یک توزیع متقارن اسست. کشیدگی یک توزیع از رابط (۱۳) به دست می‌آید. برای 
مقایسۀ دو جمعیت باید از ضریبی استفاده شود که بعد نداشته باشد مثلا" برای مقايسة دو جمعیت نمی‌توان از انحراف 
معیار استفاده نمود. زیرا ممکن است دو صفت از لحاظ واحد اندازه‌گیری یکی نباشد. برای این منظور از 


مشسخص کننده‌های نسبی پراکندگی استفاده می‌کنيم. مشخص‌کننده‌های نسبی پراکندگی از نسبت یک مشخص ES‏ 


1 Variance 

2 Coeficient of skewness 
3 Coeficient of kurtosis 
4 Standard deviation 


Lal ae ۱۷۳‏ و مخاطرات محبطی شمارة اول 


پراکندگی به یک مشخص کنندۂ مرکزی هم بعد و یا از نسبت یک مشخص SES‏ پراکندگی به یک مشخص BS‏ 
پراکندگی دیگر هم بعد به دست می‌آید. از نظر کارایی این چهار روش برآوردکردن را می‌توان به‌صورت صعودی به 
ترتیب گرافیکی» کمترین مربعات روش گشتاورها و روش بیشینۀ درستنمایی ردیف نمود. با وجود بیشترین کارایی 
استفاده از روش بیشینۀ درستنمایی در عمل بسیار مشکل است. این چهار روش با کارایی متفاوت. پارامترهای یک 
توزیع را از یک Bye‏ خاص داده‌ها محاسبه می‌کنند. این نمونه ممکن است شاخصی از جامعة اصلی باشد یا نباشد و 
استفادة از پارامترهای برآورد شده ممکن است باعث اریب (کج شدگی) نتایج شود. روش جدید دیگری که در سال 
۰ بوسیلۀ دکتر هاسکینگ از مرکز تحقیقات 18.11 ابداع شد روش گشتاورهای خحطی است که ذیلاً" در این مورد 
بحث می‌شود. 

۵-۲-۲-۲ روش گشتاورهای L‏ 

گشتاورهای 1 و گشتاورهای وزنی احتمالی شبیه به گشتاورهای عادی هستند که به‌صورت خلاصه توزیع‌های 
تتوریکی احتمالی برای نمونه‌های مشاهداتی تعریف می‌شود. مشابه با ایجاد گشستاورهای عادی» گشتاورهای 1 نیز 
می‌توانند برای اندازه‌گیری پارامترهاه فواصل برآوردها و تست مجانبی قابل استفاده باشسند. اگرچه تثوری و کاربرد 
گشتاورهای L‏ موازی با گشتاورهای رایج اسست. اما گشستاورهای L‏ دارای چندین اختلاف مهم نیز می‌باشند 
اندازه‌گیرهای نمونة گشتاورهای L‏ هميشه ترکیبات خطی از مشاهدات تصادفی هستند. آن‌ها کمیت‌هایی هستند که 
نسبت به تولید گشتاورهای gale‏ کمتر اریب هستند. این به خاطر آن اسست که اندازه‌گیرهای ایجاد کنندۀ گشستاور 
عادی به‌مانند 52 و G‏ بایستی به ترتیب مربع و مکعب مشاهدات باشند که سبب می‌شود به آن‌ها وزن بیشتری داده 
شود همراه با مشاهداتی که دور از ميانگین هستند که باعث ارب بودن واقعی و واریانس می‌شود. گشتاورهای L‏ 
می‌توانند گشتاورهای وزنی - احتمالی تعریف شوند. 

L PWMs‏ گشتاورهای وزنی احتمالی را می‌توان به‌صورت زیر نمایش داد: 

B, =Elx[F,(0] | O£) 

که FX)‏ تابع توزیع تجمعی (dÀ)‏ برای × اسست. گشستاورهای وزنی - احتمالی مورد انتظار x‏ مطابق با توان‌های 
F(x)‏ هستند. وقتی 0= ۲ باشد By‏ میانگین جریان یا :۸/2 است و اندازه‌گیرهای گشتاور L‏ در ساده‌ترین تعریف توابع 
حطی از اندازه‌گیرهای PWMS‏ هستند. اندازهگیرهای Byes‏ غیر اریب 1725 برای هر توزیعی می‌تواند از روابط 


1 n 
by =m=— > Xo 9 
11 j=l 


11-1 n -j 


سال یازدهم بررسی مقایسه‌ای بهترین روش برآورد بارش حداکثر .... 


(VY) 
al (n— j\n—j—-1\n-j-2) 
b, = 2۳ (A) 
eee 
کوچک‌ترین مشساهده اسست‎ Xi که بزرگ‌ترین مشساهده و‎ XC) جریان‌های منظم رودخانه‌ای را با‎ XG) که‎ 
کلی زیر نمایش داد:‎ doles نمایش می‌دهد. اندازه‌گیرهای 7 را می‌توان به‌صورت‎ 
n=- j 
E r 
۳ ET kL, 


(V4) 
2 


برای هر توزیعی» چهار گشتاور ابتدائی ا می‌تواند به‌آسانی با استفاده از PWMS‏ به دست آید: 
A, = Bo‏ 


(۳۰( 
 < 2, -b (۳‏ 2 
رثا + ,6 - p:‏ 6حر A‏ 
A <30 B; - 308, +12 - 6 (YY)‏ 
انداژه گیرهای نمونة خهار کشنتاور Aa yD sal‏ اندازه‌گیزهاق تمونه 8 از معادله (V4)‏ به دنست می ایند و کر 
معادلات 28 Ly yt‏ که در معادلات (۲۰) تا (۲۳) نمایش داده شده است قرار می کیرد معادلات (۲۰) تا (۲۳) 


(YY) 


به‌صورت Alles‏ عمومی زیر نمایش داده می‌شود: 
۳ 7 
(rt)‏ زمره 


مشابه با ایجاد نسبت‌های گشتاور L‏ ضریب تغییرات == ۰٤۷‏ چولگی yy‏ کشیدگی » Hosking‏ )+194( 
u‏ 


نسبت‌های گشتاور LL‏ تشریح کرد: 
(L-Cv) (Yo)‏ ضریب تغییرات nat =L‏ 
(L- Skewness) (۳۹0‏ چولگی To -7t =L‏ 
(L- Kurtosis) (v)‏ کشیدگی T, Hep‏ 
2 


٩ 3 =‏ به ترتیب ضسریب تغییرات (LCV) L‏ 


که 1,....,4 = ٤ Ay‏ گشتاور ابتدایی بآهستند و 
چولگی بآو کشیدگی ا هستند. اولین گشتاور L‏ معادل با میانگین جریان است ( *) از این نسبت‌های بی‌بعد 


Lal ae ۱۷‏ و مخاطرات محبطی شمارة اول 


می‌توان برای مقایسۀ دو جمعیت استفاده کرد. لازم به ذکر است که در این تحقیق برای برآورد پارامترهای توابع توزیع 
آماری از روش‌های گشتاور معمولی و پيشینة درست نمایی استفاده شد و نتایج استخراج شده از آن با استفاده از 
نرم‌افزار Hyfa‏ با روش گشتاور L‏ مقایسه شد. سوالی که به‌طور طبیعی پیش می‌آید این است که برای داده‌های یک 
نمون؛ خاص کدام‌یک از توزیع‌ها را به کار ببریم برای انتخاب بهترین توزیع فراوانی برای داده‌های موجود از 
آزمون‌های مختلف مانند آزمون مربع کای یا کای اسکوثر, آزمون کلموگراف-اسمیرنوف» آزمون حداقل مربعات 
استفاده می‌شود. از مان این روش‌هاء روش حداقل مربعات يا مجموع مربعات باقیمانده که در این پژوهش 
مورداستفاده قرار گرفته است را می‌توان با استفاده از فرمول زیر به دست آورد. آن توزیعی مناسب‌تر است که دارای 


مجموع le ye‏ باقیماندة کوچک‌تری باشد. 


274 ۲ 
(Qn -2,) (YA)‏ > 
R.S .S=| =‏ 
n-m‏ 
در این فرمول داریم: 
R.S.S‏ مجموع Cle‏ باقیمانده 0 مقدار برآورد شده متغیر 
,0 مقدار مشاهده شدهٌ متغیر n‏ تعداد نمونه 


به‌منظور استفاده از این روش لازم است ;ع یا دبی بر آورد شده از روش گشتاور L‏ برای هریک از ایستگاه‌ها 
وجود دارد استفاده می‌شود. به دلیل محاسبات طولانی این روش برای هریک از توزیع‌های آماری استفاده شده در این 
تحقیق» برنامه کامپیوتری در نرم‌افزار EXCEL‏ نوشته شد که این برنامه‌ها درنهایت مجموع مربعات باقیمانده برای 
هریک از ایستگاه‌ها را محاسبه می‌کنند. 


۳- نتایج و بحث 

جهت دستیابی به اهداف این پژوهش در ابتدا داده‌های ایستگاه‌های انتخاب شده هراز در دو مشسترک آماری 
برای انجام تجزیه‌وتحلیل آماده شدند. با توجه به اینکه هدف اصلی این پژوهش مقایسۀ روش‌های برآورد پارامترهای 
آماری جهت محاسبه بارش حداکثر سالیانه از لحاظ کارایی آن‌ها و داشتن حداقل خطا در محاسبه آن‌ها بوده است لذا 
در ابتدا از نرم‌افزار Hyfa‏ به‌منظور آنالیز فرکانس بارندگی‌های سالیانة ایستگاه‌های منتخب استفاده شد. در این 
پژوهش به دلیل اينکه توزیع نرمال معمولاً تطابق خوبی با روش گشتاورهای L‏ دارد موردمحاسبه قرار نگرفت 
ولی در نرم‌افزار 13910 و به روش گشتاورهای عادی و بیشینۂ درستنمائی محاسبه شد. این مراحل برای ۳ ایستگاه 


انتخاب شده بارندگی حداکثر سالیانه انجام شد. نتایج آزمون حداقل مربعات به روش گشتاورهای معمولی و بيشينة 


سال یازدهم بررسی مقایسه‌ای بهترین روش برآورد بارش حداکثر .... 


درستنمائی در این نرم‌افزار برای مقایسه با روش گشتاورهای L‏ استخراج شده‌اند. پس از این مراحل به‌منظور برآورد 
گشتاورهای del, L‏ کامپیوتری در محیط کوئیک بیسیک نوشته شد. نحوهُ محاسبه در این برنامه بدین صورت است 
که نخست داده‌های هر ایستگاه را به‌صورت نزولی مرتب می کند و سپس نسبت‌های گشتاور بآ شامل 9 و 01 و b2‏ 
و 93 (اندازه‌گیرهای نمونة مرتبه‌های اول و دوم و سوم و چهارم (Lg tS‏ را محاسبه می‌کند و Lym‏ که همان 
تعداد داده‌های هر ایستگاه است در نظر می‌گیرد و سپس گشتاور مراحل اول و دوم و سوم و چهارم ( 
2۱« 2 + 2 ,2۸ ) که به ترتیب همان میانگین» انحراف ler‏ چولگی و کشیدگی است را محاسبه کرده و 
پس از این مرحله و 7 , و 7 Ty‏ (به ترتیب ضریب تغییرات بل چولگی 1و کشیدگی (L‏ را محاسبه می‌کند که 


این مراحل برای cla tuk‏ حداکثر AUS Le‏ ایستگاه‌ها و لگاریتم بارش‌های مذکور انجام شده است. جدول (۲) 


Lubes‏ نسبت‌های گشتاور خطی برای ایستگاه‌های منتخب موردیژوهش را نشان می‌دهد. 


حدول -Y‏ محاسبۀ نسبت‌های گشتاور خطی برای بارش‌های حداکثر GULL‏ ایستگاه‌های منتخب موردپژوهش 


پارامترهای آماری 
bO‏ 
nj‏ 
bl‏ 
b2‏ 
b3‏ 
A\‏ 
AY‏ 
ar‏ 
AL‏ 


12 
13 
T4 


ا به E‏ ا 


محاسبات بر روی آن‌ها است لذا جدول (۳) پارامترهای آماری داده‌های لگاریتمی را نشان می‌دهد. 


پنجاب-هراز 
۲۰۳۸ 
tt‏ 

۱۳/۰۹ 
AVAN 
wé 

۲۳ 
۳/۶۱ 
ENR 
iN 
YMA 
۳ 


AY 


نام ایستگاه 


محمودآباد-هراز 
DAD‏ 
tt‏ 

o) 
۱۳۹۱/۶ 
Yee 

DAD 
VAT 
Vt 


۹/۹ 


کرسنگ_-هراز 
4+£/0A‏ 
tt‏ 
0Y‏ 
10/01 
۲۸۷/۳۹ 
۹۹/۸ 
۱۳/۸۹ 
Y\/08‏ 
YVOA‏ 
AY‏ 
“NA‏ 


۰/۳۲ 


مترهای آماری توزیع‌های لگاریتمی نیاز به به دست آوردن لگاریتم داده‌ها و 


جغرافیا و مخاطرات محیطی 


Ailes 


جدول ۳- محاسبة ثسبت‌های گشتاور خطی برای لگاریتم بارش‌های حداکثر GUL‏ ایستگاه‌های منتخب 


موردپژوهش 
نام ایستگاه 
پارامترهای آماری پنجاب-هراز محمودآباد-هراز کرسنگ_-هراز 

۳/۳۹۵ ۳/۰۰ ۳/۹ bO 
t£ 33 31 nj 

۱/۵۰ ۱/0۳ ۱/۸ bl 
۱/۱ ۱/۳۲ ۰/۸۰ b2 
۰۸۳ “WwW ۰/۱ 93 
۲/40 ۳/۰۰ ۲/۹ ۱ 
۰/۹ i A ۲ 
eY ۹ =Y ۳ 
۰/١ ۰۵ g A٤ 
۲۰ Nt YÉ T2 
I) ۳ ۰ 3 
۰۱۹4۸ ۰-۱۲ AA t4 


در مرحلة بعد به‌منظور انتخاب بهترین توزیع فراوانی از روش گشتاور WL‏ روش مجموع مربعات باقیمانده یا 
5 استفاده شد... پس از محاسبة RSS‏ برای AIS‏ بارش‌های ساليانة ایستگاههای متخب نتایج آن برای مقایسه با 
روش‌های گشتاورهای Spare‏ و بیشینۀ درستنمائی استخراج شدند. شکل(۲) مقایسة بین کمترین مجموع مربعات 
باقیمانده برای بارش حداکثر سالانة ایستگاه‌های منتخب به وسیل ۵ توزیع مورداستفاده در این تحقیق را با سه روش 
گشتاور معمولی. بيشينة درستنمائی و گشتاور خطی را نشان می‌دهد. با توجه به اینکه آن توزیع و روشی در این 
کل رین Ue‏ که ROSS Ch‏ ستو لا pera‏ شود ک کر کل روش os‏ کار تون 
و Hy‏ درستنمائی در ایستگاه‌های منتخب برازش بهتر و مجموع مربعات خطای کمتری را دارا است البته Job‏ 
یادآوری شود که در بعضی از توزیع‌های آماری مانند پیرسون نوع سوم در ایستگاه محمودآباد- هراز روش بيشينة 
درستنمائی Lleol‏ برازشی را با داده‌ها نشان نمی‌دهد ولی روش لوگ پیرسون نوع سوم با همین روش و در همین 


ایستگاه دارای Cues‏ مربعات Vii‏ انیت که دارای ily‏ خوبی با داده‌ها rants‏ 


سال یازدهم بررسی مقایسه‌ای بهترین روش برآورد بارش حداکثر .... wy‏ 


4 
3 
3 
3 
J 

eS 


نده 


توزیع و روش 


محمودآباد-هراز Mila Kan SM‏ پنجاب-هراز 8 
شکل -Y‏ مقایسه کمترین مجموع مربعات باقیمانده برای بارش حداکثر سالیانه ایستگاه‌های منتخب رودخانه هراز 


-rf ۶ ss ۰ 8 «€‏ . 
به سه روش گشتاور معمولی» بیشینه درستنمائی و گشتاور خطی 


در مرحلۀ بعد با استفاده از نتایج این نمودار و با توجه به کمترین مقدار RSS‏ توزیع و روش مناسب برای هر 
ایستگاه به‌صورت شکل (۳) آورده شده است. 


Lal ae \VA‏ و مخاطرات محبطی شمارة اول 


4 11.94 12.31 

N “ 6.11 
1 g 

بیشینه پیشینه درستنمائی‌6 | گشتاور معمولی ج 
درستنمائی ۱۳3 LN3‏ 7 
۳ 

محمودآباد-هراز ا کرسنگ-هراز پنجاب-هراز 3 
3 


بهترین توزیع و روش ایستگاههای منتخب هراز 


شکل ۳- بهترین توزیع و روش انتخاب شده برای بارش حداکثر سالانه از روی کمترین مجموع مربعات باقیمانده 
(R.S.S)‏ برای ایستگاه‌های منتخب رودخانه هراز 


۶- جمع‌بندی 

همان‌گونه که از شکل (۳) مشخص می‌گردد در ایستگاه پنجاب- هراز که در قسمت بالادست رودخانه و در 
ارتفاعات کوه‌های البرز قرار دارد توزیع لوگ نرمال سه متغیره و روش گشتاور معمولی با مجموع مربعات باقیمانده 
۶ بهترین روش و توزیع آماری بوده اسست. با توجه به اینکه در این بخش از رودخانه رژیم بارشی غالب» برف 
است؛ لذا این مسئله قابل‌توجه است. یکی از دلایلی که روش گشتاور خطی در این مناطق و دیگر مناطق تطابق 
مناسبی را با داده‌های بارش حداکثر سالیانه نشان نمی‌دهد. احتمالاًمی‌تواند به خاطر نوسانات بارشی سالیان اخیر و 
انیت a‏ قاده‌های elas Coley ile chy:‏ :فده است بارش atl CG Sus,‏ 
متخب موردپژوهش دارای تغییرات زیادی باشد و نتیجه آن اریب زیاد پارامترها و تأثیر روی آن‌ها مانند میانگین, 
انحراف معیار و چولگی داده‌ها بوده است. در ایستگاه کرسنگ- هراز که تقریباً در نتهای رشته‌کوه البرز از تهران به 
سمت مازندران و تقریباً در ابتدای دهانۀ جلگة مازندران در شسهرستان آمل قرار دارد نیز توزیع گمبل و روش بيشينة 
درستنمائی با عدد ۱۲/۳۱ توزیع و روش غالب این منطقه و کمترین مجموع مربعات خطا بوده است. با توجه به رژیم 
بارشی غالب این منطقه که به‌صورت بارندگی است این موضوع قابل تحقیق و پژوش است. در ایستگاه محمودآباد- 
هراز که در بخش انتهائی Hayy‏ هراز و در محل اتصال این رودخانه به دربای خزر قرار دارد نیز از بررسی ۵ توزیع 
مورداستفاده در این پژوهش شامل توزیع‌های لوک نرمال دو co pte‏ لوگ نرمال سه متغیره» پیرسون نوع سوم لوگ 
پیرسون نوع سوم و گمبل با استفاده از سه روش گشتاور معمولی» بيشينة درستنمائی و گشتاور حطی. توزیع لوگ 
پیرسون نوع سوم و روش بیشینۀ درستنمائی با کمترین مجموع مربعات خطا یعنی ۷۱۱ توزیع و روش غالب بوده 
است. با توجه به قرارگیری این ایستگاه در منطقه‌ای کاملاً جلگه‌ای که کشت غالب آن برنج است و رژیم بارشی اکثر 


سال یازدهم بررسی مقایسه‌ای بهترین روش برآورد بارش حداکثر .... ۷۹ 


مواقع سال آن بارندگی است و با توجه به تغییرات اقلیمی و نوسانات بارشی سال‌های اخیر این موضوع قابل‌بررسی 
و پژوهش است. به‌منظور مقایسۀ بهتر توزیع‌ها از روش نمره دهی به بهترین روش و توزیع نیز استفاده شده است و 
در این مقایسه به توزیعی که کمترین مجموع مربعات باقیمانده را دارد نمرۀ ۱ و به توزیعی که بیش‌ترین مجموع 
مربعات باقیمانده را دارد نمرژ 6 داده شد. همچنین به توزیعی که اصلاً برازش خوبی با داده‌ها نداشته است نمرةٌ صفر 
داده شده است و مجموع این نمرات برای هریک از توزیع‌ها به دست آمد. به‌جز نمرات صفر که مربوط به روش و 
توزیع‌های بدون تطابق ات در eb‏ موارد آن توزیعی که کمترین نمره را دارد بهترین توزیع منتخب است که نتایج 


آن به‌صورت JRE‏ (4) ارائه شده است. 


توزیع و روش 


محمودآباد-هراز ۷ کرسنگ-هراز 8 پنجاب-هراز E‏ 


شکل ۶- رتبه‌بندی نمرات کمترین مجموع مربعات باقیمانده برای بارش حداکثر سالیانه ایستگاه‌های منتخب 


رودخانه هراز به سه روش گشتاور معمولی. بیشینه درستنمائی و گشتاور خطی 


از بررسی این شکل مشخص می‌شود که در ایستگاه پنجاب- هراز در مورد روش گشتاور معمولی توزیع پیرسون 
نوع سوم و بعدازآن توزیع لوگ نرمال سه متغیره دارای برازش بهتری با داده‌های بارش سالیانه می‌باشند و در مورد 
روش بیشینۀ درستنمائی نیز توزیع پیرسون نوع سوم و پس‌ازآن توزیع گمبل در رتبه‌های بعدی برازش مناسب قرار 
دارند و در مورد روش گشتاور خطی نیز کمترین مجموع مربعات خطا مربوط به توزیع پیرسون نوع سوم و پس‌ازآن 
توزیع لوگ پیرسون نوع سوم قرار دارد. در ایستگاه کرسنگ- هراز در مورد روش گشتاور معمولی توزیع لوگ 


پیرسون نوع سوم و بعدازآن توزیع لوگ نرمال سه متغیره GLP‏ برازش بهتری با داده‌های بارش سالیانه می‌باشند و در 


Lal ae \As‏ و مخاطرات محبطی شمارة اول 


مورد روش بیشینۀ درستنمائی نیز توزیع گمبل و پس‌ازآن توزیع پیرسون نوع سوم برازش مناسب‌تری با داده‌ها را 
نشان دادند و در مورد روش گشتاور خطی کمترین مجموع مربعات خطا مربوط به توزیع پیرسون نوع سسوم و 
پسازآن توزیع لوگ پیرسون نوع سوم قرار دارد. در ایستگاه محمودآباد- هراز در مورد روش گشتاور معمولی توزیع 
پیرسون نوع سوم و بعدازآن توزیع گمبل برازش بهتری را با داده‌های بارش سالیانه نشان دادند و در مورد روش 
پيشينة درستنماتی نیز توزیع لوگ پیرسون نوع سوم گمبل و پس‌ازان توزیع لوگ نرمال سه متغیره برازش بهتری را با 
داده‌ها را نشان دادند و در مورد روش گشتاور خطی نیز کمترین مجموع مربعات خطا مربوط به توزیع پیرسون نوع 
سوم و پس‌ازآن توزیع لوگ پیرسون نوع سوم قرار دارد. با توجه به نتایج بالا کارایی هریک از روش‌های محاسبة 
پارامترهای آماری و هریک از توزیع‌ها در ایستگاه‌های منتخب منطقه مشخص شد و معلوم نمود که برای استفاده از 
آن‌ها باید دقت بیشتری نمود تا روش و توزیعی را انتخاب نمود که دارای خطای کمتر و برازش بهتری با داده‌های 
بارش سالیانه اسست. به‌منظور دست يافتن به بهترین روش و توزیع از میان نتایج UL‏ برای هریک از ایسستگاه‌های 
منتخب موردپژوهش در ilo pa‏ بعد مجموع نمرات برای هریک از روش‌ها و توزیع‌ها به دست آمد که نتایج آن در 


شکل(۵) آمده است. 


محمودآباد-هراز کرسنگ-هراز 
نام ایستگاه و روش مورد استفاده 


گشتاور خطی # بیشینه درستنمائی الا گشتاور معمولی 8 


شکل ۵- مجموع نمرات کمترین مجموع مربعات باقیمانده برای بارش حداکثر سالیانه ایستگاه‌های منتخب 


+ ۰ ۰ و ۹ e arf s oe ra‏ 
رودخانه هراز به سه روش گشتاور معمولی. بيشینه درستنمائی و گشتاور خطی 


همان‌گونه که از نتایج شکل )0( مشخص است در ایستگاه پنجاب- هراز در قسمت کوهستانی حوزة آبخیز هراز 


دارای بهترین برازش با داده‌های بارش سالیانه بوده است و روش گشتاورهای معمولی و خطی با دارا بودن نمرةٌ کل 
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۵ به‌صورت مساوی در رتبه‌های بعدی قرار دارند. در ایستگاه کرسنگ- هراز در ابتدای las‏ جلگۀ مازندران در 
رودخانة هراز در بین توزیع‌های آماری مختلف روش بیشینۀ درستنمائی دارای مجموع نمرهُ ۱۰ و بهترین برازش با 
داده‌های بارش ساانه بوده است و پسازآن به ترتیب روش گشتاور معمولی با مجموع نمرة ۱۶ و روش گشتاور 
خطی با مجموع bpa‏ ۱۵ قرار دارند که این موضوع نیز BLAS‏ غالب بودن روش بیشینه درستنمائی از نظر دارا بودن 
کمترین مجموع مربعات خطا در ایستگاه مذکور است. در ایستگاه محمودآباد هم که در بخش انتهائی رودخانهٌ هراز و 
در محل اتصال این رودخانه به Glue‏ خزر قرار دارد از بررسی مجموع مربعات خطا برای توزیع‌ها و روش‌های 
موردپژوهش نیز روش Mee‏ درستنمائی با به دست آوردن مجموع نمرة ۳ بهترین روش منتخب در این ایستگاه بوده 
است و پس‌ازآن روش گشتاور معمولی و روش گشتاور خطی به ترتیب با مجموع نمرات 1 و ۱۵ قرار دارند. نتایج 
این پژوهش را می‌توان با یافته‌های ساير محفقان در دیگر مناطق دنیا مورد مقایسه قرار داد. ازجمله پیل و 
همکاران(؟ ۲۰۰) با استفاده از روش گشتاورهای عطی علل تفاوت‌های منطقه‌ای و اقلیمی» تغییرات رواناب سالانه را 
در مناطق معتدله و خشک آفریقایی جنوبی بررسی کردند و تفاوت میزان بارش درصد پوشش گیاهی و درجه 
حرارت را علل عمده این تفاوت‌ها اعلام داشتند و همان‌گونه که این پژوهش نیز در اقلیم مرطوب شمال کشور انجام 
شده است که می‌تواند تأثیر این عوامل را در انتخاب روش و توزیع مناسب منطقه‌ای تائید نماید. در پژوهشی دیگر 
جینگای و هال" (۲۰۰۶) روش گشتاور خطی را در آنالیز منطقه‌ی سیلاب حوزه آبخیز رودخانه مینگ مورداستفاده قرار 
دادند و با استفاده از این روش مناطق همگن و توزیع پیرسون نوع سوم را به‌عنوان توزیع آماری مناسب مشخص 
کردند. با توجه به نتایج پژوهش حاضر پیشنهادهای bd‏ می‌تواند برای سایر محققین راهگشا باشد: 

۱ - پیشنهاد می‌شسود که در به‌کارگیری و استفاده از روش‌های مختلف محاسبة پارامترهای آماری برای مناطق 
مختلف احتیاط‌های لازم صورت گیرد و با توجه به تغییرات اقلیمی سالیان اخیر و تأثیر آن بر روی نوسانات بارشی 
مناطق مختلف که می‌تواند روی نتایج تحقیق SE‏ شایان توجهی را داشته باشد؛ ازاین‌رو انجام پژوهش‌های بیشتری 
در این زمینه موردنیاز است. 

۲-طول دور آماری مورداستفاده در این پژوهش دور آماری مشترک ٤٤‏ ساله بوده است و با توجه به اینکه 
نوسانات بارشی ساليانة ایستگاه‌های مختلف تأثیر زیادی را در انتخاب توزیع و روش مناسب داشته است» لازم است 
در ساير مناطق با طول دور آماری متفاوت موردبررسی قرار گیرد تا اثر آن در انتخاب توزیع مناسب به‌صورت 


واضح‌تری مشخص گردد. 


1 Peel 
2 Jingyi and Hall 


Lal ae \AY‏ و مخاطرات محبطی شمارة اول 


۳-ذر pa teal‏ ای اقغاب بهترین کر آماری از روق alll Gilg pases‏ ااه he‏ ابت 
توصیه می‌گردد آزمون‌های دیگر انتخاب بهترین توزیع آماری مانند آزمون کلم وگراف-اسمیرنوف» آزمون کای 
اسکوئر نیز IK‏ گرفته شده تا تأثیر آن بر انتخاب بهترین توزیع آماری مشخص شود. 

t‏ منطقۀ موردپژوهش رودخانة هراز در بخش مرکزی کوه‌های البرز بوده است و نظر به اينکه این رودخانه در 
Ob =‏ خود از چند بخش مختلف با آب‌وهوای گوناگون کوهستانی و جلگه‌ای عبور می‌کند تا به دریای خزر برسد و 
با توجه به قرارگیری بخشی از این رودخانه در کوهستان البرز با رژیم بارشی غالب برف و بخش دیگری در 
پائین‌دست رودخانه با رژیم بارشی غالب باران بوده است که تأثیر زیادی در انتخاب توزیع و روش مناسب داشته 
است؛ پیشنهاد می‌شود این روش مقایسه‌ای در حوزه‌های آبخیز دیگر ایران با آب‌وهوای متفاوت نیز موردبررسی و 
مقایسه قرار گیرد تا تأثیر آن در انتخاب روش و توزیع آماری مناسب برای بارش‌های حداکثر سالانه ارزیابی شود. 

تشکر و قدردانی 

از سازمان آب منطقه‌ای استان مازندران به خاطر همکاری‌های صمیمانه‌شان جهت angi‏ آمار و اطلاعات منابع آبی 


مورداستفاده در این پژوهش. کمال تشکر و سپاسگزاری را دارم. 
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